Commensurate and Incommensurate Vortex Lattice Melting in Periodic Pinning 
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We examine the melting of commensurate and incommensurate vortex lattices interacting with 
square pinning arrays through the use of numerical simulations. For weak pinning strength in the 
commensurate case we observe an order-order transition from a commensurate square vortex lattice 
to a triangular floating solid phase as a function of temperature. This floating solid phase melts into 
a liquid at still higher temperature. For strong pinning there is only a single transition from the 
square pinned lattice to the liquid state. For strong pinning in the incommensurate case, we observe 
a multi-stage melting in which the interstitial vortices become mobile first, followed by the melting 
of the entire lattice, consistent with recent imaging experiments. The initial motion of vortices in the 
incommensurate phase occurs by an exchange process of interstitial vortices with vortices located 
at the pinning sites. We have also examined the vortex melting behavior for higher matching 
fields and find that a coexistence of a commensurate pinned vortex solid with an interstitial vortex 
liquid occurs while at higher temperatures the entire vortex lattice melts. For triangular arrays at 
incommensurate fields higher than the first matching field we observe that the initial vortex motion 
can occur through a novel correlated ring excitation where a number of vortices can rotate around 
a pinned vortex. We also discuss the relevance of our results to recent experiments of colloidal 
particles interacting with periodic trap arrays. 

PACS numbers: 74.60.Ge, 74.60.Jg 
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I. INTRODUCTION 

Vortices interacting with periodic pinning arrays have 
attracted growing interest with the advances of nano- 
lithographic techniques in which the size and shape of 
individual pinning sites, as well as the geometry of the 
pinning array, can be carefully controlled. These pinning 
arrays can be composed of ID periodic modulations [Q, 
micro- holes [^]-|| or magnetic dots [p 11 in which the 
holes or dots have a diameter d significantly less than 
the dot spacing a. The holes and dots suppress the or- 
der parameter, allowing the flux lines in the mixed state 
to lower their energy, and thereby pinning the vortices. 
Strong commensurability effects in these systems are ob- 
served at the matching holds, when the number of par- 
ticles equals an integer multiple of the number of div- 
ots ]^-f4 JjLGJI . Transport measurements near T c above 
the first matching field indicate that some of the vortices 
are only weakly pinned in the interstitial regions. These 
vortices become mobile at a lower applied current than 
the vortices at the pinning sites ||. Direct evidence for 
commensurate states and interstitial vortices has been 
obtained using Lorentz force microscopy 0] and scan- 



ning Hall probe microscopy [^,^|. Magnetization mea- 
surements by Baert et al. j^] interpreted the loss of pin- 
ning above the first matching field near T c as due to a 
coexistence of an interstitial vortex liquid with a pinned 
vortex solid at the pinning sites. 

In addition to the technological applications of super- 
conductors which require strong pinning, vortices inter- 
acting with periodic pinning arrays are also an ideal re- 
alization of an elastic lattice interacting with a periodic 
substrate which can be generalized to a wide variety of 
condensed matter systems including atoms adsorbed on 
a surface |l3| , colloidal crystals interacting with interfer- 
ing laser beams |14 , and vortices in Josephson-junction 
arrays Jl5| , ^6| . A very similar system to vortices in peri- 
odic arrays that has recently been experimentally realized 
is colloidal particles interacting with optical trap arrays 
]l7| p0| . Recent experiments in these systems for colloids 
interacting with a square array observed an interesting 
transition from a commensurate square colloidal crystal 
to a triangular colloidal crystal for increasing colloidal 
densities when the colloidal interactions began to domi- 
nate over the substrate [ p^| . 

Although studies of the effects of temperature have 
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been conducted on the dynamics of vortices interacting 
with periodic pinning arrays [ pl| , the static melting has 
not been investigated. Here one could expect that vortex 
lattice melting behavior may be significantly different at 
commensurate vs incommensurate fields. Melting of vor- 
tices in wire networks has been studied for systems 
for filling fractions far less than one with Monte-Carlo 
simulations. The Monte-Carlo simulations of vortices in 
wire networks at low filling factors by Franz and Teitel 
p5[ have shown a transition from a pinned solid to a float- 
ing solid phase, both of which have the same symmetry 
as the underlying pinning lattice. For higher tempera- 
tures this floating solid phase was found to melt into a 
liquid. In the wire network case there are no interstitial 
regions, and at the matching field the floating solid phase 
is absent |fl5|| . The wire network system has several dif- 
ferences from vortices in periodic arrays in that in the 
former the vortices can be considered to move in an egg- 
carton potential. For vortices in periodic arrays the vor- 
tices would move in a muffin-tin substrate in which there 
are flat potential regions between the minima. In this 
system at low temperatures at commensurate matching, 
particles in the overlying lattice will sit in the divots or 
energy minima, while at the incommensurate matching, 
some particles will sit in the interstitial regions between 
the divots fp|j6[[T^|. Even at the matching field it can be 
possible for vortices to move into the interstitial regions 
when the vortex- vortex interactions become dominant or 
thermal fluctuations are strong enough. 

Another difference between vortices in wire networks 
compared to vortices in periodic pinning arrays appears 
for fields greater than the first matching field. The addi- 
tional vortices in the imaging experiments by Harada et 
al. j5j and simulations by Reichhardt et al |l2] ] indicate 
that in periodic pinning arrays a wide variety of intersti- 
tial vortex crystals can be stabilized at different matching 
fields. 

The vortices located in the interstitial regions can have 
significantly different static, dynamic and thermal prop- 
erties than the commensurate particles, and the coupling 
between the two different species can give rise to inter- 
esting behavior. Recent imaging experiments by Grig- 
orenko et al [Q using high resolution scanning Hall probe 
microscopy near the matching field indicate that, for 
fields slightly below and above the first matching field, 
over time the vacancies and interstitials appear to show 
site to site hopping. The nature of the hopping dynam- 
ics is however not known, such as whether the interstitial 
vortices jump from one interstitial site to another or if 
the interstitials undergo an exchange process with the 
vortices at the pinning sites. These experiments were 
done for square arrays and it is not known how the vor- 
tex hopping might occur in different geometries such as 
triangular pinning arrays. 

The melting of interstitial vortices which coexist with 
strongly pinned vortices has also been proposed by Radz- 



ihovsky |22j for vortices interacting with randomly placed 
columnar defects where there are more vortices than pin- 
ning sites. 

In this work we report the study of two-dimensional 
melting of vortices interacting with a square muffin tin 
pinning potential with the use of molecular dynamics 
simulations. For the commensurate case at low substrate 
strengths we observe a novel order-order transition from 
a square pinned vortex lattice with long-range order to 
a triangular ordered floating solid phase with quasi-long 
range order. Our simulations also allow us to examine 
the vortex dynamics at the transition into the floating 
solid where the vortices show a collective ID avalanche 
motion along the pinning rows. For increased substrate 
strength the floating solid phase is lost and the pinned 
solid melts directly to a liquid. In the incommensurate 
case for vortex densities slightly above the first match- 
ing field for strong substrates we observe a multi-stage 
melting in which the initial vortex motion occurs by an 
exchange process of interstitial and pinned vortices. 

We have also examined the effects of temperature for 
higher matching and non-matching fields greater than 
one. At the higher matching fields we find that the vor- 
tices form a solid at low temperatures. For higher tem- 
peratures we observe a coexistence of an interstitial vor- 
tex liquid with a pinned vortex solid in agreement with 
the experiments of Baert et al. At high incommensurate 
fields the extra vortices or vacancies in the interstitial 
vortex liquid become mobile. We observe the same gen- 
eral behaviors for triangular pinning arrays; however, at 
the incommensurate fields the motion of extra vortices 
or vacancies can occur by novel collective ring excita- 
tions where a group of vortices rotate around a vortex in 
a pinning site. These excitations occur when an incom- 
mensurate number of vortices cage a vortex located at a 
pinning site. 

Our system should not only be of relevance for vortices 
and colloids interacting with periodic substrates but also 
for the square to hexatic vortex transitions observed in 
certain superconductors |23| - |2"6f| , where little is known 
about the dynamics at the transitions. 

II. SIMULATION 

We model a 2D system with periodic boundary con- 
ditions in x and y with N v vortices interacting with a 
square array of N p pinning sites. We use a Langevin 
simulation for vortices at finite temperature using the 
model of Brass and Jensen (27jj2^] . We consider the vor- 
tex motion to be overdamped and integrate the following 
equation of motion. 

Here fi is the total force acting on vortex i and r/ is the 
damping coefficient. The repulsive interaction of vor- 
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tex i with the other vortices is f™ = Ylj=i /o^i(l r i — 
Tj\/X)rij. Here, A is the penetration depth and K\{r/X) 
is a modified Bessel function which falls off exponentially 
for r > A, allowing us to place a cutoff in the inter- 
actions at r = 6A for computational efficiency. The 
pinning is modeled as N p parabolic traps with f" p = 
-E^i(/pA P )(|ri-ri P) J)e(r p -|r i -r^|)f^. Here, 6 is 
the Heaviside step function, r p is the range of the pinning 
potential, f p is the maximum pinning force of the wells, 

f y = (ji - Tj)/\Ti - rj\ and r$ = (r. t - r{ p) )/l r i - r^]. 
The thermal force has the properties < ff{t) >= and 
< fi{t)ff{t ) >= 2T)k B T6ij6(t-t'). In this work we set 
f} = ks = 1- All lengths, fields and forces are given in 
units of A, $o/A 2 and /o respectively. For most of the 
results shown here the density of vortices is 0.25<5>o/A 2 . 
The initial ground state of the system is found using sim- 
ulated annealing. We then slowly increase T and ana- 
lyze the coordination number, P n (from the Voronoi or 
Wigner Seitz cell construction), and the structure factor, 

5(k) = (l/L 2 )^y k -[ r *- r ^ 

where L is the system size length. 

Another measure we employ to characterize the ther- 
mally induced transition is the sum of the displacements 
of the vortices from their initial positions at T = to 
their positions at a higher T: 

N v 

d„,(T) = ^|r i (T)-r i (0)| 2 

i=l 

If the vortex lattice configuration is the same as the initial 
state than d w — 0. If a portion of the vortices move but 
then become frozen again d w will jump up. If the vortex 
lattice becomes molten d w will increase to a saturation 
value of L/2. 

To visualize vortex motion we plot the trajectories that 
the vortices follow for a period of time. To examine the 
dynamic response in the phases we can add a driving 
term representing a Lorentz force from an applied 
current to Eq. 1 and sum over the net vortex velocities 
V x = Y^jZi v i 'X. To analyze finite size scaling in S(k) we 
conduct simulations with the same vortex and pin densi- 
ties for system sizes with a length L for systems ranging 
from L = 24A to L = 60A with N v = 144 to N v = 900. 

III. COMMENSURATE MELTING 

In Fig. 1 we present the evolution of the coordination 
numbers Pq and P4, the average power of the secondary 
peaks in the structure factor < S(k) >, and the vor- 
tex displacements d w for increasing T for a system at 
the matching field with L = 48 A. Here f p = 0.03/o, 
r p = 0.17A and a = 2.0A. For < T < 0.0012 the pinned 



vortex lattice has the same square symmetry as the pin- 
ning lattice, as indicated by P4 ~ 1.0. The initial increase 
in d w for T < 0.0012 is due to the vortices moving ran- 
domly within pinning sites but remaining confined to a 
distance r p . Near T = 0.0012 there is a structural transi- 
tion in the vortex lattice from square lattice to triangular, 
indicated by the sudden drop in P4 to ~ and the large 
increase in Pq to 0.98. In Fig. 1(b) we plot the average 
< 5(k) > of the appropriate secondary peaks in S(k). 
For T < 0.0012 we observe only four secondary peaks in 
S'(k) as indicated in the left inset in Fig. 1(a). At the 
square-triangular transition there is initially a large drop 
in < 5(k) > from 0.95 to 0.2 and then a recovery to 0.55. 
The position of the secondary peaks shifts at the tran- 
sition from four-fold to six-fold order. During the tran- 
sition a mixture of both four and six- fold peaks appear, 
reducing the overall power at each peak. In the triangu- 
lar phase there is almost no power for the four-fold order. 
The square to triangular transition can also be observed 
as a jump in d w (Fig. 1(c)), coinciding with the drops 
in P4 and < S(k) >, which is due to the shifting of the 
vortex positions during the transition. In the triangular 
lattice phase, d w remains roughly constant, indicating 
that diffusion is not occurring and that the vortex lattice 
is still a solid. Near T = 0.0049 the triangular vortex 
lattice melts into a liquid, indicated by the increase in 
d w as well as the drops in < S(k) > and P§. We note 
that Pq will not drop to zero for a liquid state since even 
for random particle distributions a portion of the vortices 
will always have six- fold coordination numbers. We have 
also conducted simulations where we cool down in tem- 
perature starting at T = 0.008. In this case we do not 
observe any hysteresis in the liquid to triangular phase; 
however, hysteresis is present in the triangular to square 
transition, suggesting that the transition is first order in 
nature. We have also examined the response of the phase 
to an applied drive and find that only the square pinned 
phase has a finite critical depinning force. 

The finite size scaling behavior of the structure factor 

S(k)/L 2 ~L- e 

which is described more fully in Jl5| gives us further 
information about the nature of the phases. For long 
range order, 9 = 0, for floating solids or hexatic phases 
< 9 < 2.0, and for the liquid phase 9 = 2.0. For 
different system sizes we average the power of the sec- 
ondary peaks in < S(k) > for 150 frames. In the inset of 
Fig. 1(b) we plot < 5(k) > /L 2 (where N v ~ L 2 ) versus 
N v for the three different phases. In the commensurate 
phase (T = 0.0005) we find 9 w 0.0 corresponding to the 
pinned lattice. For T = 0.003, in the triangular phase, 
9 = 0.33 ± 0.04 which is consistent with a floating solid 
phase. In the floating solid phase we find a small varia- 
tion of 9 with T which we will address elsewhere. In the 
liquid phase (T = 0.007) we find 9 = 1.98 ± 0.04. 
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A. Dynamics of the Order-Order Transition 

In order to gain insight into the dynamics of the square 
to triangular transition, in Fig. 2 we show the real space 
images as well as the individual vortex trajectories in the 
pinned solid, the pinned-solid to floating solid transition, 
and the liquid phase for the system in Fig. 1. In the 
pinned solid phase shown in Fig. 2(a) (T = 0.0005) we 
find a square vortex lattice with little or no detectable 
vortex motion. Fig. 2(b) (T = 0.00135) shows the dy- 
namics at the beginning of the pinned to floating solid 
transition. Here alternate rows of vortices move in a sud- 
den ID avalanche where the vortices in one row move in a 
collective manner while the other rows remain static. All 
the vortices in a given row shift together as the vortices 
jump out of the pinning sites and move to the interstitial 
regions. The motionless rows shown in Fig. 2(b) also be- 
gin to move at a slightly later time. This shifting of the 
vortex positions produces the jump in d w at T — 0.0012 
seen in Fig. 1(c) at the square to floating solid transi- 
tion. In Fig. 2(c) (T = 0.0015), in the floating solid 
phase, the vortex lattice has triangular order and shows 
increased thermal wandering around the equilibrium po- 
sitions as can be seen in the vortex trajectories. This 
increased thermal wandering also smears the secondary 
peaks in S'(k), reducing the maximum value to 0.55, as in 
Fig. 1(b). In Fig. 2(d) (T = 0.0065), in the molten state 
the vortices are disordered and diffuse in a random man- 
ner corresponding to the large value of d w in Fig. 1(c). 

B. Phase Diagram for the Commensurate Case 

By conducting a series of simulations with different 
substrate strengths and using the criteria from Fig. 1 to 
determine the onset of the different phases, we construct 
the phase diagram in Fig. 3. The range of temperatures 
over which the commensurate solid appears increases lin- 
early for increasing f p . The floating solid phase only ap- 
pears for f p < 0.125/o. The line separating the floating 
solid to liquid phase is roughly constant with tempera- 
ture for lower substrate strengths at a value that corre- 
sponds to roughly the zero pinning melting temperature, 
T m = 0.0048 (as also obtained with the diffusion and 
structure factor measurements). All three phase bound- 
aries can be understood from simple considerations. The 
floating solid to liquid line is independent of f p because 
it is determined by the vortex-vortex interactions and 
not the pinning force. For f p > 0.125/o the commensu- 
rate solid melts directly into a liquid and < S(k) > has 
the behavior shown in the inset of Fig. 3. The melting 
transition is seen to shift to higher temperatures with 
increasing pinning strength, as the melting occurs when 
thermal fluctuations are able to overcome the increasing 
pinning strength. Finally the commensuratc-to-floating 



solid follows fp^T with the low f p being washed out be- 
fore the vortex- vortex interactions are washed out. The 
nature of the floating solid to liquid transition, such as 
whether a dislocation-mediated melting transition occurs 
]l3[ , is beyond the scope of this work. 

To further compare our results with those of vortices 
in wire networks |l5| ], in which the floating solid phase 
with the same symmetry as the pinned phase is observed 
as a function of commensurability, in our case the pinned 
to floating solid transition is seen as a function of sub- 
strate strength. Further, we observe a square to triangu- 
lar transition due to the fact that a portion of the vortices 
shift into pin free regions, which is not possible in wire 
networks. 

We note that in our simulation wc imposed nearly 
square boundary conditions which can create defects or 
stress in the triangular lattice. We do not find any de- 
fects in our triangular lattice and our finite size scaling 
is consistent with that in |ll| . It is also possible that the 
square pinned solid phase may be enhanced for a higher 
temperature range by the square boundary conditions; 
however, the general trends of the phase diagram should 
still hold. 

IV. INCOMMENSURATE MELTING NEAR THE 
FIRST MATCHING FIELD 

We now consider the vortex behavior above the first 
matching field when there are more vortices than pins. 
We find that for large substrate strengths, f p > 0.25/o, a 
commensurate sub-lattice forms, with the additional vor- 
tices sitting in the interstitial regions. In Fig. 4 we show 
the melting behavior for a system with N v = 1.063-ZVp and 
f p = O.6/0. In the imaging experiments of Grigorenko et 
al. |jj , the interstitial motion was observed for a field of 
B/Bj, = 1.04. Here, for T < 0.0016 both the commensu- 
rate vortices and the interstitial vortices remain pinned 
as seen in Fig. 4(a). At these low temperatures, the inter- 
stitial vortices are pinned by the potential cage created 
by vortices at the pinning sites. At T > 0.0016 vortex 
diffusion begins to occur as indicated by the increase in 
d w in Fig. 5(b). This is due to the onset of motion of the 
incommensurations as seen in the vortex trajectories in 
Fig. 4(b). The motion is not restricted to interstitial vor- 
tices. Instead, interstitial and pinned vortices exchange 
places, and over a period of time all the vortices take part 
in the motion. The exchange process can be understood 
by considering that the interstitial vortex will produce 
a force fi n on the nearby commensurate vortices which 
effectively reduces the pinning force to \f p \ — \fin\- Since 
the melting transition is a function of the pinning force 
the melting temperature will be reduced. 

In Fig. 5(a) the onset of the incommensuration mo- 
tion is also marked by a drop in < S(k) > from 0.96 
to approximately 0.7. In contrast to the sharp drop in 
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< 5(k) > seen in the commensurate case, here < 5(k) > 
still retains strong four-fold order due to the presence 
of the background lattice of pinned vortices as seen in 
Fig. 4(b). For higher temperatures < 5(k) > gradu- 
ally decreases with most order being lost for T > 0.0175, 
corresponding to the melting temperature for the com- 
mensurate case with the same substrate strength. For 
simulations where we cool down from T > 0.0175 we do 
not observe any hysteresis. We have also conducted sim- 
ulations in which we apply a constant drive of fd = 0.012, 
gradually increase T and measure the average vortex ve- 
locity < V x >. For T < 0.0016, in the pinned interstitial 
phase < V x >= 0. At the interstitial transition there is 
a jump in < V x > indicating the onset of vortex motion. 
Finally at the overall vortex melting transition, < V x > 
jumps to a value corresponding to the entire lattice flow- 
ing. 

V. MELTING FOR HIGHER COMMENSURATE 
AND INCOMMENSURATE FIELDS 

A. Commensurate Melting 

We have also investigated the melting in square and 
triangular pinning arrays for the higher matching fields 
B/B^ = 2,3,4, 5 and 8. At these fields ordered vor- 
tex crystals are formed |||l2|. In Fig. 6 we show the 
vortex positions and trajectories for the square pinning 
array with B/B^ = 4.0 for increasing temperatures. For 
T < 0.004 the vortex lattice is ordered as seen in exper- 
iments and simulations. For T > 0.004 [Fig. 6(b)] the 
interstitial vortex lattice melts; however, the vortices at 
the pinning sites remain immobile so that a coexistence 
of a solid and a liquid phase occurs. The interstitial vor- 
tices are constrained to move in a square grid and there 
is a region around the pinning sites which the interstitial 
vortices do not enter due to the repulsion of the vortices 
located at the pinning sites. Under an applied drive the 
depinning force is finite for low temperatures but goes to 
zero at the onset of the interstitial vortex liquid transi- 
tion. This result is in agreement with the interpretations 
of the experiments by Baert et al. in which the pin- 
ning force goes to zero near T c when interstitial vortices 
are present, but is still finite for fields at which vortices 
are only located at the pinning sites. For increased T, 
as in Fig. 6(c), the vortices at the pinning sites become 
depinned as well and the entire lattice is in the molten 
state. The general behavior of the melting from Fig. 6 
is also seen at the other matching fields we have investi- 
gated with transitions from the solid, to solid-liquid co- 
existence, to the liquid state. 



B. Incommensurate Melting 

In Fig. 7 we show the melting behavior for the incom- 
mensurate field B/Bcf, = 4.1. Here the onset of vortex 
diffusion occurs for a much lower temperature than for 
B/B<p — 4.0. The vortex motion occurs at defect sites 
where additional interstitial vortices are present in the 
ordered B/B^ = 4.0 interstitial vortex crystal. These 
extra interstitial vortices are less strongly pinned than 
the other interstitial vortices. The vortex motion can 
occur by the continuous motion of a single vortex over 
a certain distance, or by a pulse like motion in which a 
series of individual vortices move by one lattice constant 
as the pulse moves through. If the trajectories are drawn 
for a longer time period the interstitial vortex diffusion 
can be seen to occur through the entire lattice. As the 
temperature is increased the vortices at the pinning sites 
will become depinned. 

C. Collective Ring Excitations for Incommensurate 
Fields with Triangular Pinning 

We have found that the melting behavior for square 
and triangular pinning arrays is similar at the matching 
fields; however, at the incommensurate fields we observe 
clear differences between the two pinning geometries. In 
Fig. 8 we show the melting behavior for a system with tri- 
angular pinning for B / B^ = 3.08. At low temperatures 
the system is frozen but as T is increased motion begins 
to occur in the form of ring excitations. Here a portion 
of the vortices rotate around a vortex at a pinning site; 
however, no net diffusion of interstitial vortices is occur- 
ring. For B/Bcj, = 3.0 each vortex at a pinning site is 
surrounded by six interstitial vortices. For fields above 
or below B/B^ = 3.0 a portion of the pinned vortices 
will be surrounded by 5 or 7 vortices, and it is at these 
sites where the ring excitations occur. Due to the sym- 
metry of the pinning array the vortices at the pinning 
sites create a potential around each pinning site which 
has six minima, so that at B/B^ = 3.0 the vortex con- 
figuration is particularly stable. When there are seven or 
five vortices around a pinning sites, they are incommen- 
surate with the six mimima, and therefore the thermal 
kicks can cause the interstitial vortices to slowly rotate. 

As the applied field is moved further away from the 
matching field, the number of incommensurations in- 
creases and the number of ring excitations increases as 
well. Such ring excitations are observed in triangular 
pinning arrays for most incommensurate fields but with 
different numbers of vortices forming the ring. As the 
temperature increases, the initial diffusion of the vortices 
occurs via a hopping mechanism in which the extra vor- 
tex in the rotating ring hops out of the ring and starts 
the rotation of a different ring. The rate of this hopping 
increases with the temperature. 
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VI. SUMMARY 

In summary, we have investigated the melting of vor- 
tex lattices interacting with periodic pinning arrays for 
both the commensurate and incommensurate case. For 
the first matching field with weak substrates we have ob- 
served a transition in the vortex lattice from a square 
commensurate solid to a triangular floating solid which 
melts into a liquid at higher temperatures. This tran- 
sition occurs by ID collective motion of vortices shift- 
ing from the pinning sites into the interstitial regions. 
For strong substrates the floating solid phase is lost and 
the vortex lattice melts directly into a liquid. For fields 
slightly above the first matching field for strong disorder 
we observe a multi-stage melting where the interstitial 
vortices are highly mobile consistent with recent imaging 
experiments. The motion of the interstitial vortices oc- 
curs by the exchange of interstitial with pinned vortices 
rather than with the interstitials hopping from one in- 
terstitial site to another. As the temperature increases 
more vortices become mobile. 

We have also examined the melting behavior for match- 
ing fields greater than one. At low temperature the in- 
terstitial vortices form a crystalline state. At higher tem- 
peratures we observe a coexistence of an interstitial vor- 
tex liquid with a commensurate pinned vortex lattice. 
The interstitial vortices can diffuse in the regions between 
the pinned vortices. We also observe a depletetion zone 
around the pinned vortices which interstitial vortices do 
not enter due to the repulsion from the pinned vortex. 
At higher temperatures the entire vortex lattice melts. 
At higher incommensurate fields the additional defects 
(extra interstitials or vacancies) in the interstitial vor- 
tex crystal become mobile and diffuse in the interstitial 
regions. With triangular pinning arrays at incommen- 
surate fields the initial vortex motion occurs through a 
novel collective ring excitation where interstitial vortices 
can rotate around a pinned vortex. These rotations occur 
when an incommensurate number of interstitial vortices 
surround a pinned vortex. 

We briefly discuss some experimental systems in which 
these effects could be observed. For vortex lattices in 
periodic pinning arrays, vortex lattice melting is most 
relevant to high Tc superconductors. Our results are 
only for a 2D system whereas 3D effects can be rele- 
vant to melting. Our model would thus best describe 
high-temperature superconductors with periodic colum- 
nar defects where the vortices would have line like be- 
havior. However, several of the results here should still 
be relevant for low temperature superconductors. For the 
low temperature superconductors melting would only oc- 
cur very near T c . For strong pinning such as arrays of 
holes, a floating solid phase at the first matching field 
would not be observable and the system should go di- 
rectly from the pinned solid to the normal phase with 



a very small region of a disordered vortex lattice. For 
weak pinning such as small magnetic dots or weak defect 
arrays the pinned solid to floating solid transition should 
occur well below T c and should be observable for the low 
temperature superconductors. For a low temperature su- 
perconductor the phase diagram in Fig. 3 would be mod- 
ified, with the floating solid to liquid transition replaced 
by the floating solid to normal transition, although again 
a small vortex liquid state may occur just below T c . Also 
the pinned solid to vortex liquid transition would remain 
flat rather than increasing with f p . The solid to float- 
ing solid transition could be imaged with scanning Hall 
probe measurements, neutron scattering, or Lorentz mi- 
croscopy. In addition transport measurements would also 
be able to reveal the loss of pinning at the transition. 

The vortex behavior at the incommensurate fields 
should also be visible in the low-temperature supercon- 
ductors. The recent imaging experiments of Grigorenko 
et al. 0] have already found evidence for the motion of 
the highly mobile interstitial vortices in these types of 
systems. These same imaging techniques should be able 
to image the pinned vortex lattice coexisting with the 
interstitial vortex liquid. In triangular pinning arrays 
collective ring excitations would appear as a smeared 
ring around pinning sites with the Hall probe arrays and 
the dynamics could be directly imaged with Lorentz mi- 
croscopy. 

Our work is also directly relevant to colloidal particles 
interacting with periodic pinning. A particularly promis- 
ing realization of this system is optical trap arrays in 
which the pinning strength can be easily tuned. In such 
a system our phase diagram can be directly tested. Re- 
cent experiments have already seen evidence for a square 
to triangular transition as the colloid density is increased 
fl9f . In addition the dynamics of the interstitial colloids 
can be directly imaged with video microscopy to deter- 
mine if the interstitials hop directly from one interstitial 
site to another, or if they move through an exchange pro- 
cess with a pinned colloid. In addition it should be also 
be possible to observe the melting behavior for higher 
matching colloidal densities and the collective ring ex- 
citations for colloids interacting with triangular optical 
trap arrays. 
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FIG. 1. (a) The fraction of four-fold P4 (thick line) and 
six-fold P§ coordinated (thin line) vortices versus T. (b) The 
average power of the secondary peaks, < S'(k) > /N v versus 
T. The position of the peaks shifts from four-fold to six-fold 
near T = 0.001. The peaks are lost at T = 0.005. The inset 
of (b) shows the scaling of the secondary peaks for different 
system sizes where N v ~ L 2 . For the pinned phase (circles), 
T = 0.0005 and S(k) ~ L°; for the triangular lattice (trian- 
gles), T = 0.003 and S(k) ~ L~ 0,33 ; and for the liquid phase 
(squares), T = 0.007 and S(k) ~ L" 1 ' 98 . (c) The displace- 
ments d w /N v =< |r(T) — r(0)| 2 > of the vortices versus T. 



FIG. 2. The vortex positions and trajectories for a 
24A x 24A section of a 48A x 48A system for (a) the pinned 
solid (T = 0.0005), (b) the beginning of the pinned solid to 
floating solid transition (T = 0.00135), (c) the floating sold 
phase (T = 0.0015), and (d) the liquid phase (T = 0.0065). 
In (b) alternate rows of vortices shift positions. For slightly 
higher temperatures the other vortices also shift in position. 
This shifting is reflected in the jump in d w at T — 0.0012 in 
Fig. 1(c). 




0.002 0.004 0.006 0.008 



FIG. 3. The phase diagram of f p versus T. The open 
circles represent the pinned solid to floating solid transition. 
The filled circles represent the transition into the liquid phase. 
The inset shows the behavior of < S(k) > /N v versus T for 
f p = 0.125/o where the pinned solid melts directly to a liquid. 
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FIG. 4. (a) A snapshot of a 10A x 10A region of the sample 
with B/Bcf, — 1.06. The vortex positions and trajectories for 
T = 0.0005 show that some vortices are pinned at interstitial 
positions, (b) At T — 0.0036 the motion of incommensura- 
tions can be seen while the background lattice remains pinned. 
The vortex trajectories in (b) show that vortex motion occurs 
by the interstitial vortices pushing other vortices off the pin- 
ning sites which in turn move into interstitial regions and 
repeat the process. 
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FIG. 5. The behavior of < S(k) > /N v and d w /N v versus 
T for a system with f p = 0.6/o and N v = 1.063JV P . Motion 
in the vortex lattice starts at T = 0.0016 as indicated in the 
increase in d (b). 
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FIG. 6. (a) The vortex positions (black circles) and pin- 
ning sites (open circles) for B/B<f, = 4.0 for a system with a 
square pinning array at T = 0.001. Here the vortex system 
is in a crystalline state, (b) For T = 0.004 the vortex trajec- 
tories show that the interstitial vortices are in a liquid state 
while the vortices at the pinning sites remain immobile, (c) 
At T — 0.01 all the vortices are mobile. 
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FIG. 7. The vortex positions (black circles) and pinning 
sites (open circles) for a system with the same pinning param- 
eters as in Fig. 6 for B/B^ = 4.1. Here T = 0.001 and the 
vortex trajectories show a portion of the interstitial vortices 
are mobile. 



* 



® * * ® J ® £ C ' 

f © © ST© © © 
. vl) © V © © 




* 


© • 






* 












© 


© 



FIG. 8. The vortex positions (black circles) and pinning 
sites (open circles) for a system with same parameters as in 
Fig. 6 with a triangular pinning array with B /B$ — 3.08 and 
T — 0.001. Here the vortex motion can be seen to occur 
in a ring excitation where vortices can rotate around pinned 
vortices. The ring excitations occur where there are seven 
vortices around a pinned vortex. A second type of excitation 
can be seen in the trianglular trajectories where three vortices 
can rotate. 
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